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NON-HOLONOMIC CONSTRAINED SYSTEMS 
AS IMPLICIT DIFFERENTIAL EQUATIONS 



Abstract. Non-holonomic constraints, both in the Lagragian and Hamiltonian formalism, 
are discussed from the geometrical viewpoint of implicit differential equations. A precise 
statement of both problems is presented remarking the similarities and differences with 
other classical problems with constraints. In our discussion, apart from a constraint sub- 
manifold, a field of permitted directions and a system of reaction forces are given, the later 
being in principle unrelated to the constraint submanifold. An implicit differential equa- 
tion is associated to a non-holonomic problem using the Tulczyjew's geometrical description 
of the Legendre transformation. The integrable part of this implicit differential equation 
is extracted using an adapted version of the integrability algorithm. Moreover, sufficient 
conditions are found that guarantees the compatibility of the non-holonomic problem, i.e., 
that assures that the integrability algorithm stops at first step, and moreover it implies the 
existence of a vector field whose integral curves are the solutions to the problem. In addition 
this vector field turns out to be a second order differential equation. These compatibility 
conditions are shown to include as particular cases many others obtained previously by other 
authors. Several examples and further lines of development of the subject are also discussed. 



1 Introduction 

Non-holonomic constraints have been the subject of deep analysis (not exent from 
some controversy) since the dawn of Analitical Mechanics. In fact D'Alembert's 
principle of virtual work jnj, p5| , and Gauss principle of least constraint jl7| can 
be considered to be the first solutions to the analysis of systems with constraints, 
holonomic or not. A golden age for the subject came with the contributions of O. 
Holder ||, G. Hamel @, P. Appell §, E. Delassus @, T. Levi-Civita @, N. 
G. Chetaev ||, etc., when the discussion of (linear) non-holonomic constraints in 
Lagrangian mechanics was considered sistematically. After the quantum revolution 
this classical problem was kept frozen in the limbo of the "postponed" problems. 
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However, there has been papers that from time to time have addressed some of the 
weak points of the discussions as they were left in the thirties, for instance, the 
existence and uniqueness of solutions, the inadequacy of Chetaev's conditions, etc. 
Some contributions in this transition era can be found in the papers by R.J. Eden p"5[ , 
V. Valcovici p3fl , J. Neimark and N. Fufaev p5[ ], R. van Dooren ]T3], V. Rumiantsev 
Y. Pironneau §6|, etc - 

Almost simultaneously with this period a quite revolution was taking place re- 
garding the foundations of the old discipline of Mechanics. Geometry was used in a 
sistematic way to set up its foundations and new and old ideas from Geometry and 
Mechanics were fusioning in an harmonic and pleasant picture that has been evolving 
until today and is not yet completely finished. W. Tulczyjew was one of the 
pionnccrs of this Geometric revolution and his ideas and insight have inspired many 
developments in the field as is reflected for instance by the variety of contributions 
to these Proceedings. Sooner or later the geometrical shock wave had to reach also 
the problem of non-holonomic constraints. In fact, the first geometrical description of 
non-holonomic constraints took place in relatively early times in a paper by L. D. Fad- 
deev and A. M. Vershik |Q. The importance of this paper was not recognized until 
quite recently even if it contained the first general result on the existence of dynamics 
of Lagrangian systems with (not necessarily linear) non-holonomic constraints. Even 
previous to this, J. Klein had already addressed the problem of constraints using the 
geometry of Lagrangian systems p3[ . Much recent are the papers by E. Massa et al 
discussing the geometrical meaning of Chetaev's conditions |3^]; G. Giachetta and 
M. de Leon et al [^8) using jet bundle geometry techniques; J. Koiller J24|, L. Bates et 
al [^| and Bloch et al |3| discussing non-holonomic constraints with symmetry. More 
contributions can be found in the papers by L. Cushmann et al Q] using a symplec- 
tic splitting to describe the dynamics in presence of non-holonomic constraints; J. 
Carihena & M. F. Ranada || geometrizing Lagrange's multipliers; R. Weber Gf|, A. 
Van der Shaft et al Q, P. Dazord |I) and C. M.-Marle Q addressing the problem 
from the Hamitonian viewpoint; Sarlet et al |4C| ] using connection theory; de Leon 
et al using almost product structures and projectors, F. Barone et al Q using 
Tulczyjew's ideas to set the frame for generalized Lagrangian systems, F. Cardin et 
al combining the vakonomic and non-holonomic approach, and finally CM. -Marie 
in these Proceedings revisiting Faddeev- Vershik conditions. 

Our modest contribution to this old subject will consist in using Tulczyjew's idea 
of modelling mechanical systems as implicit differential equations j30), J33|, to 
discuss non-holonomic constraints. This approach was already taken by S. Bencnti 
H who was able to obtain a set of sufficient conditions for the existence of dynamics 
for linear non-holonomic constraints and apply it successfully to solve the problem 
of constrained geodesies. In our approach we will compare first the problem of La- 
grangian systems with non-holonomic constraints with the problem of constrained 
Hamiltonian systems. We will set up the geometry of Lagrangian and Hamiltonian 
systems with non-holonomic constraints and after a brief discussion of implicit differ- 
ential equations and the integrability algorithm, we will apply it to the non-holonomic 
problem. We will find that a small modification of the integrability algorithm is 
needed to encompass the restrictions imposed by non-holonomic constraints to the 
solutions of the problem. From the analysis of the adapted integrability algorithm 
some particular cases arise immediately. The simplest non-trivial case is discussed 
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thoroughly and a general condition for the existence of a solution, that turns out 
to be a vector field, is discussed. This geometrical condition contains as particular 
cases Faddeev-Vershik's condition, Chetaev's conditions, Bates' regularity condition, 
Benenti's conditions, Marie's conditions, etc. 

We will apply these ideas to discuss a several examples that are inspired in old 
models like Appell's machine, etc. 

2 Non-holonomic constrained systems 

The discussion to follow of non-holonomic constraints will be set in the realm of 
Tulczyjew's triple @, 

T*(TQ) T(T*Q) T*(T*Q). (2.1) 

We will recall that a is the Tulczyjew's canonical symplectomorphism from T(T*Q) 
(with its canonical symplectic structure ujq) to T*(TQ) (with its canonical structure 
cotq), and (3 is the canonical symplectomorphism defined by the former symplectic 
structure. 

The problem of constrained Hamiltonian sytems consists in determining the equa- 
tions of motion for a system specified by a constraint submanifold C C T*Q and 
a Hamiltonian function H.C — *■ M. The submanifold C is usually, but not always, 
determined by a self-consistency condition of the system under study and it is often 
the result of a Dirac's type constraint algorithm |jl3), fisj] . The submanifold 

D = {veT c {T*Q)\uj Q (v,u) = (dH,u), VueTC}, (2.2) 

defines an implicit differential equation on T(T*Q), called a Dirac system, whose 
analysis provide the solution of the posed problem. Such analysis was succesfully 
done in |m| and a set of necessary and sufficient conditions for its integrability was 
given in the following theorem. 

Theorem 1 The Dirac system D defined by a constrained Hamiltonian system is 
integrable iff C is coisotropic and H projects along the characteristic distribution of 
C. 



Despite the similarities between the Hamiltonian constrained problem and the 
problem of non-holonomic constraints, the latter has a different nature. In the non- 
holonomic Lagrangian problem a Lagrangian L will be given in all TQ that defines 
the unconstrained or "free" system. In fact, we can assume in what follows that L 
is regular. This means that on TQ we have a dynamical vector field given by the 
Euler-Lagrange vector field T l defined by the equation, 

C TL Q L = dL, (2.3) 

or equivantly {L regular), 

ir L uJL = dE L , 
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where ®l — S*(dL) is the Poincare-Cartan 1-form defined by L, ujl = — d&L is the 
Cartan 2-form and El — A(i) — L is the energy of L, with A the Liouville vector 
field on TQ. 

The canonical tensor field S on TQ is defined in local natural coordinates (q 1 , v" 1 ) 
on TQ by 

S = ^- <g> dq\ 

and it allows us to characterize SODE's as S(T) = A. The kernel and image of S 
consists of vertical vector fields. On the other hand S acts by duality on forms and 
the kernel and image of S* consists on horizontal 1-forms. A remarkable property of 
S and ujl is given by the formula, 

i s uJL = 0, 

or cquivalently, 

S* o ui L = ° S, (2.4) 
where Cjl denotes the map T(TQ) — > T*(TQ) defined by contraction with wj,. 

Notice that L being regular implies that the 2-form ujl is nondegenerate. Conse- 
quently we can define its inverse bivector as, 

K L {a,l3)=u L {ul x {a),u7 1 {P)l Va,/3 G T*(TQ). (2.5) 



Now for reasons that in principle have no relation with the Lagrangian L, a sub- 
manifold C C TQ will be selected, the constraint submanifold of the problem. If the 
submanifold C is not of the form TN with N C Q, or more generally, if C is not the 
total space of an integrable distribution T> defined along a submanifold N c Q, we 
will say that the constraints are non-holonomic. In general will not be tangent to 
C, this is, the Lagrangian system defined by L will evolve in time without keeping 
within the limits imposed by C. If we want the system to remain in C then Tl must 
be changed. For that we will assume that there is a set of "forces" F C T*(TQ) 
that allow us to act upon the dynamical system Tl and eventually to make it to be 
confined to C. 

The modified systems that we can obtain from by means of the forces F are 
given by the family of vector fields, 

rer L + A L (F), (2.6) 



where is the Poisson tensor defined by the Cartan 2-form u>l, eq. (2.5). Notice 



that Euler-Lagrange's equations (2.3) in the presence of external forces / = fi(q, v)dq l 
are modified as 

£r® L =dL + f, 



thus, spanning the Lie derivative, we will get that if / G F, then T verifies eq. (2.6). 

If we want the resulting system to be a system of "mechanical" type, the system 
of "forces" will have to be given by horizontal 1-forms, i.e., S*(F) = 0. Then a simple 
computation shows that T is a SODE because S(Al(F)) = Al(S*(F)) = where we 



have used eqs. (2.4) and (|2.5|) 
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Figure 1: Non-holonomic constraints data, field of directions and forces. 



The system of forces F will be obtained cither by a detailed analysis of the con- 
straint submanifold or they will be given in an independent way. For instance, a 
large class of non-holonomic constraints are originated by the interaction between the 
surfaces of different components of the system. In this category fall sliding, rolling 
and friction constraints, which are supposed to create linear relations between the 
velocities of the components of the system. However, we can also imagine that our 
system is subjected to the action of servomechanisms or other devices, that modify its 
dynamical state in such a way that certain "a priori" given conditions are satisfied, 
for instance, limitations in the acceleration of the center of mass of the system. The 
system of forces F can be postulated from the non-holonomic constraints, i.e., they 
will be supposed to be created by the constraints and they will be supposed to have an 
explicit relation with them. This is the case of the so called Chetaev's conditions for 
non-holonomic constraints (see Section 5). Some of these possibilities will be explored 
later on. In this sense the system of forces will be called reaction or control forces, 
depending if they are derived from the constraints or are imposed externally. 

Thus the Lagrangian non-holonomic constraint problem can be stated as follows: 
given a regular Lagrangian L, a constraint submanifold C and a system of forces F, 
determine wether or not there exists a vector field T of the form given by eq. (2.6) 
parallel to C. 



This problem admits a simple generalization that is of interest in a variety of 
situations. We can replace the symplectic manifold (TQ,ljl) by a Poisson manifold 
(P, A) and the dynamical vector field Tl by a Hamiltonian vector field Th = A(dH). 
We will consider a constraint submanifold C C P and, contrary to the situation 
discussed above, we can think that a further restriction on the possible directions of 
the vector field Th can be imposed, i.e., a field of allowed or permitted directions D 
along C will be introduced^. The geometrical model for it will consists in a vector 
subbundle p: D — > C of rp\TP — > P (see Fig. [|). The reaction forces now will 
be modelled by a subbundle r\: F — > C of np:T*P — > P. Then, the generalized 
non-holonomic problem can be stated as follows: 



Poisson non-holonomic constraint problem: Given a Hamiltonian vector field 
Th on a Poisson manifold P with Poisson tensor A, a constraint submanifold C, a 
field of permitted directions D and a system of control forces F, determine if there is 
a vector field T of the form 

rer ff + A(F), (2.7) 

1 If we are discussing a Lagrangian system, such restriction can be understood as a limitation on 
the accelerations of the system. 
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contained in D. 

We should point it out that in both versions of the non-holonomic constraint 
problem the sought vector field can exists or not, it can exists only on a subset of 
the constraint submanifold, and even if it exists, it is not necessarily unique. It is 
clear that all kind of possibilities can actually occur and it is not difficult to exhibit 
examples of them [[33) . The discussion in Sections 3 and 4 address the existence and 
uniqueness of such dynamical vector fields. Particular situations will be also discussed 
where explicit conditions can be done that will allow to give definite answers to them. 

3 Implicit differential equations 

The main tool to analyze the general problems discussed above will be the geometrical 
setting developed in |33| to describe implicit differential equations. 

An implicit differential equation on the manifold P is a submanifold E C TP. A 
solution of E is any curve 7: 1 — ► P, I C M, such that the tangent curve (j{t), ^(t j) € 
E for all tel. The implicit differential equation will be said to be integrable at a 
point if there exists a solution 7 of E such that the tangent curve passes through it. 
The implicit differential equation will be said to be integrable if it is integrable at all 
its points. Integrability does not imply uniqueness. The integrable part of E is the 
subset of all integrable points of E. The integrability problem consists in identifying 
such subset. 

Denoting as before the canonical projection TP — > P by rp, a sufficient condition 
for the integrability of E is 

E C TC, (3.8) 

where C = Tp(E), provided that the projection rp restricted to E is a submersion 
onto C. 

Extracting the integrable part of E: A recursive algorithm was presented in [^3| 
that allows to extract the integrable part of an implicit differential equation E. We 
shall define the submanifolds 

E Q = E, C Q = C, (3.9) 

and recursively for every k > 1, 

E k = E k ^ n TCfc-i, C k = T P {E k ), (3.10) 

then, eventually the recursive construction will stabilize in the sense that E k = 
E k+ \ = ■ ■ ■ = Eoo, and C k = C k +\ = ■ ■ ■ = Coo. It is clear by construction that 
Eoo C TCoq. Then, provided that the adequate regularity conditions are satisfied 
during the application of the algorithm, the implicit differential equation Eoo will be 
integrable and it will solve the integrability problem. 

In what follows we will refer to the non-holonomic problem defined by the La- 
grangian L, the constraint submanifold C, the distribution D and the system of 
forces F, as the system (L, C, D, F). 
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4 Integrability of non-holonomic systems 



4.1 The implicit differential equation associated to a system 
with non-holonomic constraints 

The implicit differential equation associated to a non-holonomic system can be de- 
scribed in the Lagrangian and/or the Hamiltonian formalism. 

Lagrangian picture: We will be given the data defining a non-holonomic con- 
strained Lagrangian system, i.e., a regular Lagrangian L, a field of permitted direc- 
tions D C Tc{TQ) along the constraint submanifold C C TQ and a system of forces 
F C Tq(TQ). Then, the implicit differential equation associated to L,D,F will be 
defined by 

E = a- 1 {dL + F) CT(T*Q). (4.11) 

Notice that tt*q{E) — C* is the constraint submanifold considered in the problem 
and it coincides with TTl (C) where Tl denotes the Legendre transformation defined 
by L (see below). It is also noticeable that dL + F defines an affine subbundle of 
T*(TQ) along C. The regularity of L guarantees the transversality of the submanifold 
E with respect to the projection map tt*q. 



Hamiltonian picture: As in the Lagrangian case, a non-holonomic constrained 
Hamiltonian system is defined by a Hamiltonian function H on T*Q, a permitted 
directions field D C T C ,{T*Q) along the submanifold C* C T*Q and a system 
of forces F C T£,(T*Q). The implicit differential equation associated to the data 
(H,C*,D,F) is defined by, 

E — f3~ 1 (dH + F), (4.12) 

that projects onto C* . 



4.2 Non-holonomic integrability algorithm 



If we apply the integrability algorithm given by eqs. (3.9)-( 3T0|) to the implicit 



differential equation associated to a non-holonomic Lagrangian system given by eq. 
(4.11)^1 we will obtain, 



E a =E = a~ 1 (dL + F), C = t t *q(E)=C*, 
£i=£ nTC , C 1 = tt*q{Ex),..., 

and eventually c TCoq. However we must notice that the stability of the algo- 
rithm does not imply automatically that E^ is going to define a vector field lying in 
D and the integrable equation will not be in general a solution of the non-holonomic 
problem. 

Then, the integrability algorithm needs a small adaptation to render the required 
outcome in the context of non-holonomic systems. We need to translate the field of 
allowed directions to T*Q. The natural way to do that is by means of the Legendre 

2 N ntirp that a similar algorithm can be applied to non-holonomic Hamiltonian problems, eq. 
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Figure 2: The non-holonomic integrability algorithm. 



transformation T L defined by L. The derivative of this map TI L : T(TQ) -> T(T*Q) 
allows to define the bundle of permitted directions on T*Q as 

D = TT L {D) CT(T*Q), 

then we will construct an adapted integrability algorithm as follows, 

Eq = E, C — C*, d = d, 

E l = E r\b , c 1 = t t *q{e 1 ), bi = D \ Cl mrci, 
E k = e^ n D fc _i, c k = T T , Q (E k ), b k = bk-i \c„ nTCfe, 



and then, eventually when the algorithm stops, we will obtain the stable data E^ , 
Doa , Coo defining an integrable implicit differential equation (see Fig. ^). 

Notice that Eoo C -Doo C TCoo which together with the appropriate regularity 
assumptions guarantees the integrability of Eoo and its compatibility with the non- 
holonomic constraints of the problem. 

Obviously all kind of situations for Eoo, boo, Coo can happen. The simplest 
nontrivial^ situation occurs when E\ = Eoo, C* = Coo and Doo — D. Then, the 
algorithm stops in the first step and, 

E l = Enb = Er\TT L {D), Ci = T T . Q (Ei) = c*, £>! = bnrc* = tt l {d). 

(4-13) 

We will say that the regular Lagrangian L and the non-holonomic constraints D, 



F are compatible if the conditions above, eq. (4.13), are met, i.e., the adapted 
integrability algorithm for them stops at the first step. In this case, the integrable 
implicit differential equation defined by them are the submanifold of T(T*Q) given 
by, 

E 00 = a- 1 (dL + F)nTF L {D), (4.14) 



3 Nontrivial here means that C is actually a submanifold of TQ of codimension at least one and 
F is a subbundle of rank at least one. 
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and it gives the solution of the non-holonomic Lagrangian problem (L,C, D, F). 

We must notice here that the outcome of the integrability algorithm does not 
necessarily provide us with a SODE. If we require that the solution must be a SODE, 
then we must modify the algorithm again to incorporate this fact. This problem has 
been exhaustively discussed in the setting of Lagrangians systems were it is known as 
the SODE problem for singular Lagrangians. 



4.3 Compatibility conditions 

It is important to remark again that the previous characterization of compatibility, eq. 



(4.13), does not imply the uniqueness of solutions passing through its points. This is a 
general feature of non-holonomic constraints often encoutered in practical discussions. 
It is also important to notice again that even if the solutions were unique, i.e., the 
implicit differential equation Eoo would define a bona fide differential equation, they 
would not be necessarily solutions of a SODE. The following theorem gives sufficient 
conditions for the non-holonomic Lagrangian problem (L,C, D, F) to be compatible 
and, as a bonus, it is found that they also guarantee the uniqueness and the second 
order character of it. 



Theorem 2 Let (L,C, D, F) be a non-holonomic Lagrangian problem. If 

rank(F) = corank(D), (4-15) 

and 

D°nF ± =0, (4.16) 

then the Lagrangian and the non-holonomic constraints are compatible. Moreover, 

there exists a unique SODE T such that its graph is the integrable implicit differential 
equation a~ 1 (dL + F) n TTl{D) associated to it. 



Proof: The symbol D° denotes the annihilator of D, i.e., 

D° = {aeT*(TQ)\(a,D)=Q}, 

and F denotes the symplectic orthogonal to F with respect to the Poisson tensor 
Al, that is, 

F ± = {feT£(TQ) | A L (/,F)=0}. 

The vector fields constructed from to solve the Lagrangian non-holonomic problem 
have the form 

T = T L +A L (f), feF. 
If r € D, then pairing the previous equation with elements in D°, we will get 

A L (a,f) = {T L ,a), VaeD . 

Consider the map tf>: F — ► (D )* defined by <p(f)(a) = A L (a, /), for all a £ D°, whose 
kernel is F n (D )^. Notice that (F fl (D )^) 1 - = F x + D°. But dim^- 1 + D°) = 
dimF^ + dimD because of eq. ( 4.16 ), then, 

dim(F n (D )- 1 -) =2n- (dimF 1 + dimD ) = 2n - (2n - rankF) - corankD = 0. 
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Then we conclude that F H (D ) 1 - = and thus conditions ( fhl5| )-( p~l~(f ) imply the 



injectivity of the map <j). Moreover, again because of rank(F) = corank(D), the 
previous map is surjective. Hence the map <f> is a bundle isomorphism and there 
exists a unique element / G F such that <fi(f) = (Tl, .). Then there exists a unique 
vector field that gives a solution of the non-holonomic Lagrangian problem, i.e., the 
submanifold E cuts TTl(D) along a submanifold which is the graph of the Legendrc 
transfom of the vector field T. Moreover, because the forces / a horizontal 1-forms, 
the vector field A^(/) is vertical and T is a SODE. □ 

Remark: It follows from the proof of the previous theorem that the conditions for 
compatibility can also be written as 

rank^) = corank(L>°), F n (D ) 1 - = 0, 

or they can be simultaneously encoded as 

Tc(TQ) = F 1 - © D°. 

4.4 Local expressions 

Because of its practical interest we shall write the conditions in Thm. || in local 
coordinates. 

We will suppose that the submanifold C is defined locally by the set of functions 

4> a {q,v) = 0, a=l,...,r. 

Then dimC = 2n-r. The annihilator of the subbundle D C TC will verify (TC)° C 
D°. Thus, we will decompose D° — (TC)°®K. We will suppose that K is a subbundle 
with a local basis in the chart q l , v l given by the set of 1-forms 

(3 b = ^dq 1 + P\dv\ b = l,...,S. 

A system of local 1-forms generating (TC) is given by d<p a , a = 1, . . . ,r. Thus, 
D° will be locally generated by the set of 1-forms d<fi a ,(3 b , that will be collectively 
denoted by 

I3 C = ^dq 1 + (ifdv\ C = 1, . . . ,r + s. 

The corank of D will then be (r + s). 

We will now assume that the system of forces F will be defined by a subbundle 



of rank r + s, satisfying thus eq. (4.15), and it will have a system of local generating 
horizontal 1-forms f c — ffdq 1 , C = 1, . . . ,r + s. The local expression of the Poisson 
tensor A^ is given by 

a „wi d d , „.,• d d 
Al = W lJ — A — ^ + M rj — A — — r . 

oq 1 ov l av % ov^ 



The intersection condition (4.16) is equivalent to the statement that Al defines a 
nondegenerate pairing between F and D , i.e., that the matrix 

M BG = W ij (q,v)ff(q,v)pf(q,v), 
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is nondegenerate. In the particular case that D = TC, i.e., there are no limitation to 
the permitted directions for the vector field, then the set of 1-forms spanning D° are 
given by d(j) a , a = 1, . . . ,r. Then, the compatibility condition becomes, 

0^debM ab = W ii (q,v)ft(q,v)& 



Linear non-holonomic constraints: 

The situation that usually is dealt with is when the submanifold C is an affine 
subbundle A of TQ, i.e., locally, C is defined by local functions, 

4> a (q, v) = tf{q)v l + = 0, a = 1, . . . , r. 

The tangent submanifold TA C Ta(TQ) is defined by the set of 1-forms d(f> a , i.e., 
(TA)° = span{ d<fi a \ a = 1, . . . , r }. Hence, because 

namely, 




= tf(q), 



and the compatibility condition becomes, 

detM Qb = det W tJ (q, v)f?(q, v)^(q) ^ 0. (4.17) 

If, in adition we accept Chetaev's forces (see next Section), i.e, we define, 

F = S*{TA)°, (4.18) 

the subbundle F will be generated by the 1-forms /i a = n1(q)dq % because S*(d<fi a ) = 
fif(q)dq t , and the compatibility condition will be simply given by 

detM Qb = detW ij (q,v)^(q)fj, b j (q) + 0, (4.19) 

reproducing results in ]27|| . Notice that if L is a mechanical Lagrangian L = T — V , 
where T is the kinetic energy corresponding to a Riemannian metric g 1 then, 

AI ab = gVtfiq^q), 

which is obviously invertiblc. 



5 Comparison with previous results 

We have already seen in the previous section, that particular choices of the set of 
reaction forces F, for instance those given in eq. ( 4.18| ), give special expressions for 



the compatibility conditions discussed in Section 4 before, Thm. |2|. A special choice 
for the reaction forces is given by Chetaev's conditions. In intrinsic terms the system 
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of forces corresponding to a non-holonomic constraint submanifold C C TQ under 
Chetaev's conditions is given by the subbundle (see also p2|), 

F = S*((TC)°). (5.20) 

More generally, if we consider a bundle of permitted directions D — > C, then 
Chetaev's bundle of reaction forces will be given by (see also [|tJ)j 

F = S*(D°). (5.21) 

Following Jl6| we will say that D is admissible if S* is injective when restricted 
to D°, i.e., D° does not contains horizontal 1-forms. This implies that the rank of 
the Chetaev bundle coincides with the rank of D°, i.e., rankF 1 = corankD. Hence, 



Faddeev's admissibility condition is a particular instance of the condition (4.15) in 
Thm. |. 

The Hessian of the Lagrangian L can be defined as the symmetric (0, 2)-tensor^|, 

H L (a,/3)=A L (S*(a),0), a,/3eT*(TQ). 

Then, we will say that L is definite if the Hessian Hl is definite as a symmetric tensor. 
Then, if L is definite, Hl is nondegenerate when restricted to the subspace D°, but 
H L {D°,D°) = Al(S*(D°),D°) = A L (F,D°), and A L is nondegenerate on the pair F, 
D°, i.e., F 1 - n D° = 0. Then if L is definite ("normal" according to the terminology 
in |nj) this implies condition ( 4.1 6j ) . Then Thm. || implies the main result in p6[ , 



Theorem 3 Let L be a definite Lagrangian and a non-holonomic constraint defined 
by the permitted field of directions D — > C '. If D is admissible and the reaction forces 
are given by the Chetaev's bundle F = S*(D°), then there exists a unique SODE T 
solution of the non-holonomic Lagragian problem defined by L and D — > C. 



A different condition is used in |3| and |27j to find a solution of the Lagrangian 
non-holonomic problem. Now we will assume that a submanifold C C TQ is given 
and D = TC. Let Sc be the distribution 

S c = S{TC^). (5.22) 
A Lagrangian system L with a non-holonomic constraint is said to be regular if 

S c nTC = 0. (5.23) 
The following theorem is again a particular instance of Thm. ^| 

Theorem 4 If (L, C) is a regular non-holonomic Lagrangian system, then it has a 
unique solution. In addition this solution is a SODE. 

4 Notice that is symmetric because of eq. (E^l) 
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Proof: The analysis of the regularity condition Sc H TC = leads to, 



= lu l (S c n TC) = Co L (Sc) n u L {TC) = u L {S{TC^)) n Qi L {{{TCff), (5.24) 

where as before <2>l is the natural bundle map T(TQ) — > T*(TQ) defined by con- 
traction with the 2-form cj£. Then it is obvious that u>l(TC ) = TC°. In fact, 
a € ^(rC- 1 ) iff there is u € TC^ such that u L (u) = a, but u> L (u,TC) = 0, then 
a G TC and conversely. Then, continuing with the computation in cq. ( 5.24 ) we 
get, 

U! L (S(TC x ))nu! L (((TC)°) Q ) = S*(w L (TC ,x ))n(TC' ) ± = S^TC^niTC ) 1 - = Fn(D°) x , 
with F = S*{TC°) and D = TC, and again we are in the conditions of Thm. |. □ 

As it was pointed before, a similar approach to this paper was adopted by S. 
Benenti in || . There, the problem of linear non-holonomic Lagrangian or Hamiltonian 
systems was addressed and the compatibility conditions exhibited there (Props. 2 and 



3) are equivalent to the regularity of L and eq. (4.19). The nonlinear case is also 
discussed and the conditions in Prop. 4 of |f| are equivalent respectively to the 
regularity of the Lagrangian, the rank condition in Thm. and the regularity of 
the matrix M ab in eq. (4.17) plus Chetaev's conditions; hence the conclusions there 



follow from Thm. |2j. 

To end this section we will comment on the characterization given by CM. Marie 
in j3lj. The subbundles W (called the projection bundle) and TD, the tangent bundle 
of the Hamiltonian constraint submanifold, are W — TTl(Sc) and D = C*, using 
the notation in this paper. Then, W (&TD is the translation to T*Q of the regularity 



condition eq. (5.23). Then, the existence and uniqueness of solutions stated in [Bl| 



(Prop. 2.15 and Thm. 2.16) are again particular instances of the previous results. 



6 Examples and applications 
6.1 Disk rolling on a surface 

As a simple example we will discuss first the a disk rolling vertically on a rough 
surface. The disk of mass M and radius R will be described by the coordinates x, y 
of the contact point with the surface, the angle 9 defined by the plane containing 
the disk and a fixed plane normal to the surface, and the angle <f) parametrizing the 
position of the disk with respect to its center (see also Fig. ||). The configuration 
space of the system will be Q = K 2 xT 2 . The Lagrangian describing the free system 
is, 

L = l -M(x 2 + f) + l -hP + (6.25) 

whith Ii , I 2 the corresponding moments of inertia. The constraint submanifold C C 
TQ is given by the rolling conditions, 

*i = x - Rcos6(j) = 0; * 2 = y - Rsin6<j) = 0. (6.26) 
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In the cotangent bundle T*Q = T*R 2 xT*T 2 we introduce coordinates (x, y, 9, 4>,p x ,p y ,P9,Pcf>), 
and finally in T(T*Q) we will have coordinates [x, y, 9, 4>, p x ,p y ,pe,P</>; x, y, 9, <j>, p x ,p y ,P9,p<p)- 
The annihilator of TC on Tg(TQ) will be spanned by d^i, <N> 2 - Using Chetaev's 
forces, we will have, F = S*(TC)°, 

F = {X 1 (dx~Rcos9d(j)) + X 2 {dy-Rsm9d(j)) \ X U X 2 el}. 

Thus, the affine subbundle dL + F C T^(TQ) is given by, 

dL + F = { M xdx + Mydy + h<pd<j> + I 2 9d9 + 

+\i(dx - Rcos9d(j)) + X 2 (dy - Rsm9d4>) | X X X 2 G K }. 

The map a: T(T*Q) -> T*(TQ) is defined in coordinates (q l ,Pi; q\p\) for T(T*Q) 
and (q l , v l ; n,Si) by 

s i oa=p l , no a = pi, v l oa = q l , q l oa = q\ 

or equivalently, a(q l ,Pi,q l ,Pi) — (q l , q l ,pi,pi). Then, the submanifold E = a~ 1 (dL + 
F) is given by the equations, 

E = {(x,y,9,(j>,p x ,p y ,p e ,p^x,y,9,<p,p x ,p y ,pe,p4,) eT(T*Q) \ 

Px = Xi,p y = X 2 ,p (t> = -Aii?cos6» - X 2 Rsm9,pe =Q,p x = Mx,p y = My, 
p^ = h4>,pe = I 2 9,^ 1 = o,* 2 = o}. 
The Legendre transformation Fl- TQ — > T*Q is given by 

p x = Mx,p y = My^p^ = h<j>,p g = I 2 9, 
and maps the submanifold C into the submanifold C* <ZT*Q given by 

C* = { (x, y, 9, <\>;p x ,p y ,p e ,p <i> ) € T*Q \ hp x = mRcos9p ( j,, hp y = mRsin9p lj> }. 

It is clear that tt*q{E) = C* , and thus EdTC* — E\ and the integrability algorithm 
stopsQ. Computing the intersection of TC* with E we obtain, 

MR sin 9 x MR cos 9 
Ai = P4,pe, X 2 = — p^pe, p e = 0, p 4> = 0. 

In this form we will obtain the following vector field T defined on C* , 

i?cos(9 d Rsin9 d p^ d p e d MRsin9 d MRcos9 d 
h ox h v dy h dep I 2 89 I X I 2 dp x hh v dp y 

which is mapped by F^ 1 into a SODE, or equivalently the following set of equations 
of motion, 

z = ^f%, y = ^f%, * = 9=%, 

whose solutions are found immediately. 

The results thus obtained are in full agreement with the solutions obtained in any 
elementary course in Analytical Mechanics and illustrate the predictions of Thm. ||. 

5 This was known in advance because L is definite and the constraints are admissible. 
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Figure 3: Appell's machine with p = 0. 



6.2 Appell's machine 

We will discuss now the well-known Appell's machine (see Figure ||). 
p = 0: It consists in a disk of radius R with a small cylinder rigidily attached to it 
of radius r with total mass M . A rope passes through the small cylinder and hangs 
a mass m on the other extreme of a vertical frame. The frame is of negligible mass 
and the disk rolls over the surface. The coordinates of the contact point of the main 
wheel will be (x, y) and the coordinates of the small mass m will be thereby (x, y, z). 
The angle defined by the horizontal axis of the wheel with the Ox axis will be denoted 
by 9 and the angular position of the wheel will be denoted by <f> as in the example 
before. The configuration space being Q — M 3 x T 2 and the tangent bundle becomes 
the 10-dimensional manifold TQ = TR 3 x TT 2 . The La grangian of the system will 
be given by 

L = -M{x 2 + f) + -m(x 2 + y 2 + z 2 ) + h^ 2 + h 2 9 2 + mgz. 

The rolling conditions between the disk and the floor and the nonsliding condition 
between the rope and the cylinder define the constraint submanifold C C TQ. The 
7-dimensional submanifold C is then characterized by the equations, 

*i = x - Rcos6<j> = 0, * 2 = y - i? sin 00 = 0, * 3 = r <j) - z = 0. 

The non-holonomic constraints are linear and the system can be considered to be 
of Chetaev's type with F = S*(TC)°. The Legendre transformation of C gives the 
submanifold C* , 

C* = {(x,y,z,6,(j),p x ,py,p z ,pg,p$) €T*Q \ Iip x = (M + m)Rcos6p < f > , 
hp y = (M + TO)i?sin6»p0, Ixp z = mRp^ }. 

Repeating the computations of the previous section we will arrive to the following 
equations of motion on the submanifold C* C T*Q, 
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Px = Ai, p y = A 2 , p z = A 3 - mgr, (6.27) 
P0 = -X 1 Rcos9 - X 2 Rsmd - A 3 r, p e = 0. (6.28) 

Thus, solving the Lagrange's multipliers, Ai, A2, A3, i.e., computing E n TC* , we get, 
" mgr 



h + (M + m)R 2 + ■ 



= 0, 



which is integrated immediately. Two comments are pertinent here. Notice that the 
last equation is simply the conservation of the z component of the angular momentum, 
which is obvious from the analysis of the forces acting upon the system. Secondly, 
we know in advance that the system is compatible and has a unique and well defined 
solution without having to consider it as a limit case p — ► of a "well posed" system. 
In fact we will see immediately that the case p 7^ (see Fig. |J) has extra difficulties 
because of the nonconservation of angular momentum. 

p / 0: The system now consists as before on a disk an a small cylinder rigidily 
attached to it. A rope passes through the small cylinder and hangs a mass m on 
the other extreme of a frame of length p. The frame is of negligible mass and slides 
without friction. Now the coordinates of the contact point of the main wheel will be 
{xdtUd) and the coordinates of the mass m will be (x,y,z). The angle defined by 
the horizontal axis of the wheel with the Ox axis will be denoted again by 9 and the 
angular position of the wheel will be denoted by 4>. The Lagrangian of the system 
will be given now by 

L = \M (x 2 d + y 2 D ) + im(i 2 + y 2 + z 2 ) + l -I^ 2 + l -I 2 Q 2 - mgz. 

The geometry of the system imposes that x — xq = pcosO and y — yn = psin^, 
(that can be understood as holonomic constraints). We can eliminate them from the 
Lagrangian L and we get a system on Q = R 3 x T 2 as in the situation before. The 
constraints now are given by, 

*i = x D - Rcos9(f> = 0, * 2 = i)D - -Rsin6»(/> = 0, * 3 = rcj) - z = 0. 

Repeating the computations and using the Legendre transform to write now the 
equations of motion directly on C C TQ we will obtain, 

(M + m)x + M p6 2 cos 6 + M p6 sin = A x , 
(M + m)y + Mp() 2 sin 9 - Mp6 cos = A 2 , 
m'z + mg = —A3 

Ii<f> = —XiRcos9 — X 2 Rsm9 + A 3 r, 

(I 2 - Mp 2 )9 - MpR04> = \ lP sm9 - \ 2 pcos9 

Eliminating the Lagrange's multipliers, we obtain <f> and 9 variables, 

{I 2 + (M + m)p 2 )9 + pRm9<j) = 0, (6.29) 
(ii + [m + M)R 2 + mr 2 )4> - MRp9 2 = —mrg. (6.30) 
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Figure 4: Appell's machine with p ^ 0. 



Eq. (6.29) shows that there is no conservation of the angular momentum, the reason 
being that for the mass m to move rigidily with the frame, this has to act on it with a 
horizontal force. The previous equations, ( |6.30 ), ( 6.29| ), can be solved and we obtain 
<fi(t) implicitely from the integral, 



V2 



0o - mgrcf, + - + ^^ ^ z 2 -(M+^ (* ~ 0o) 



t-tn 



and 



9(t)=9 Q + I e - i2-^Z^ Wt) - M 



dt. 



la 



6.3 A variation of Benenti's problem 

Finally we will consider some variants of an example proposed by Benenti g| (see Fig. 

D- 

First we will solve the problem of two point massess forced to move on a plane 
with parallel velocities. The configuration space will be Q = I 2 x R 2 and we will 
denote by (xi, yi) the position of the particle of mass mi and by {x2, J/2) the position 
for the particle of mass mi- The constraint on the velocities is given by the function 
on the tangent bundle TQ, 

* = xifa ~ &2V1- (6.31) 

The Lagrangian of the system is, 

L = ~mi(xl + y{) + \^n 2 {x\ + yl). 

The non-holonomic constraint ^ is a genuine non-linear non-holomic constraint. If 
we solve the system with Chetaev's forces ( 5.20| ) we will obtain the following system 
of equations on TQ, 

rnxxx = Xy 2 , miy\ = -\x 2 , m 2 x 2 = -Xy%, m 2 y2 = \x 1 , 
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Figure 5: Benenti's system. 



thus derivating eq. (6.31), and substituting on it, we get, 

mi m 2 J 

Hence, A = and Chetaev's forces vanish. The solution of the system is simply 
free motion of the two particles with initial parallel velocities. Notice that Chetaev's 
forces have the form fx = X(y2dxi — x 2 dyi), fi = X(—yidx 2 + Xidy?) where fi is the 
force acting on the particle i, thus the only solution leaving the submanifold $ = 
invariant is obtained with A = 0. However if we choose the system of forces 

F = { Ximifidxx + mifsdyi + m 2 fidx 2 + m-ifidy-}) | A G R}, 

then we will have the equations, 

*i = A/i, t/i = A/ 2 , x 2 = Xfi, j/2 = A/ 2 , 

thus if the initial velocities were parallel, the corresponding motions will be a solution 
for arbitrary A. Thus, in this particular situation there are non unique solution of the 
problem. 

The problem can be modified substituting the point masses by disks rolling without 
sliding on the plane. Then we will be considering two copies of the system discussed 
in Section 6.1., eqs. ( |6.25| ), ( J6.26 ) , with the additional constraint, 

$5 = ViAv 2 , (6.32) 

where denotes the velocity of the center of mass of the it h dis k. Replacing the 
constraints given by the rolling conditions, eq. ( |6.26 ), on cq. ( 6.32 ), we will obtain, 

RnR 2 sin(0 2 - 8i)fafa = 0. 

Then, either one of the disks is still (0j = 0) and the angles 8i are arbitrary or 
if the disks are rolling ^ 0), then 6\ = 9 2 = 6 and the nonlinear constraint 
^5 is redundant. Then the two disks roll independently and freely keeping parallel 
directions. 
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If we remove the rolling conditions and we keep only the parallelism condition 
\&5 = 0, then we will obtain again the solutions of the point masses problem before, 
i.e., free motion with initial parallel velocities. Under this condition, the disks slide 
without rolling. It is obvious that using non-Chetaev forces, other solutions can be 
found. 

7 Conclusions and outlook 

We have reviewed the theory of Lagrangian systems with non-holonomic constraints 
from the viewpoint of implicit differential equations in the realm of symplectic and 
tangent bundle geometry. From the beginning we have deliberately separated the non- 
holonomic constraints from the forces or controls that we can use to force the system 
to satisfy them, emphasizing in this way that non-Chetaev's systems can be included 
in this picture. In fact, we have found general conditions that guarantee the existence 
and uniqueness of solutions for this generalized non-holonomic problem. We have 
especialized this general compatibility conditions to different particular situations 
already discussed in the literature. 

The discussion has been made in the autonomous case but it is obvious that it 
can be extended to the time-dependent setting easily (see for instance Q , p8| , Q ) . 
Some of the beautiful aspects of Tulczyjew's triple are lost in a time-dependent setting 
and this justifies to keep the discussion at the autonomous level. The usual extended 
phase space trick can be used to include time as an additional variable and proceed 
as it has been done here. An alternative path would be to extend Tulcyjew's triple to 
the setting of cosymplectic/contact geometry to recover the same geometrical setting 
without introducing spurious variables. 

Subtler is the problem of extending the theory of non-holonomic constraints to 
singular Lagrangians. This is a significant extension of the theory because physi- 
cal Lagrangians are very often singular. Singular lagrangians introduce their own 
constraints, that arise as integrability conditions of an implicit differential equation. 
A simultaneous analysis of Lagrangian constraints and non-holonomic constraints is 
needed to obtain the equations of motion of the theory. This problem has been ad- 
dressed in H and ||. 

An important aspect of the theory of non-holonomic constraints is its relation 
with Quantum Mechanics. Apart from a paper by R.J. Eden there has been lit- 
tle attention to the quantization of systems with non-holonomic constraints (see also 
[|37|). One reason for this is that fundamental theories do not include non-holonomic 
constraints because non-holonomic constraints are phenomenological models that of- 
fer a reasonable description for the (often unknown) true fundamental theory that 
describes the interaction between the surfaces of the system in contact, or other in- 
teractions present in the system which cannot be described at the same fundamental 
level. Some aspects concerning the Hamiltonian structure of non-holonomic systems 
and its quantization are being discussed in p2| . 
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